Nuclear mid-infrared properties of nearby 
low-luminosity AGN 



D Asmusi, S F Honig^, P Gandhi^, A Smette^ and W J Duschl^'S 

^ Institut fiir Theoretische Physik und Astrophysik, Christian- Albrechts-Universitat zu Kiel, 

Leibnizstr. 15, 24098 Kiel Germany 
^ Department of Physics, University of California in Santa Barbara, Broida Hall, Santa 

Barbara, CA 93106-9530, USA 

Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency, 
3-1-1 Yoshinodai, chuo-ku, Sagamihara, Kanagawa 252-5210, Japan 
^ European Southern Observatory, Casilla 19001, Santiago 19, Chile 

^ Steward Observatory, The University of Arizona, 933 N. Cherry Ave, Tucson, AZ 85721, 
USA 

E-mail: asmusSastrophysik . uni-kiel . de 

Abstract. We present ground-based high-spatial resolution mid-infrared (MIR) observations 
of 20 nearby low-luminosity AGN (LLAGN) with VLT/VISIR and the preliminary analysis of 
a new sample of 10 low-luminosity Seyferts observed with Gemini/Michelle. LLAGN are of 
great interest because these objects are the most common among active galaxies, especially in 
the nearby universe. Studying them in great detail makes it possible to investigate the AGN 
evolution over cosmic timescale. Indeed, many LLAGN likely represent the final stage of an 
AGN's lifetime. We show that even at low luminosities and accretion rates nuclear unresolved 
MIR emission is present in most objects. Compared to lower spatial resolution Spitzer /IRS 
spectra, the high-resolution MIR photometry exhibits significantly lower fluxes and different 
PAH emission feature properties in many cases. By using scaled Spitzer /IRS spectra of typical 
starburst galaxies, we show that the star formation contribution to the 12 /im emission is minor 
in the central parsecs of most LLAGN. Therefore, the observed MIR emission in the VISIR 
and Michelle data is most likely emitted by the AGN itself, which, for higher luminosity AGN, 
is interpreted as thermal emission from a dusty torus. Furthermore, the 12 /im emission of the 
LLAGN is strongly correlated with the absorption corrected 2-10 keV luminosity and the MIR- 
X-ray correlation found previously for AGN is extended to a range from 10*° to lO^'^erg/s. 
This correlation is independent of the object type, and in particular the low-luminosity Seyferts 
observed with Michelle fall exactly on the power-law fit valid for brighter AGN. In addition, no 
dependency of the MIR-X-ray ratio on the accretion rate is found. These results are consistent 
with the unification model being applicable even in the probed low-luminosity regime. 



1. Introduction 

One of the major discoveries of the last two decades in extragalactic astronomy was that the 
structure and evolution of the central kiloparsec region of galactic nuclei are strongly connected 
to the properties and evolution of the central super-massive black hole (SMBH) that presumably 
resides in the center of each galaxy. In particular, the growth of the SMBH seems to be related 
to major galaxy mergers leading to various phases of nuclear activity (see contribution by Duschl 
et al.). These active galactic nuclei (AGN) can be described by the unification model of AGN [T] 



that contains the SMBH in the nucleus, surrounded by an accretion disk and, further outward, 
a dusty obscuring torus- hke structure. This torus presumably blocks the direct view onto the 
central region depending on the viewing angle of the observer, therefore explaining the different 
AGN sub-classes observed. However, for a more accurate description, various modifications to 
this standard model may arguably be necessary (e.g., as suggested in the contributions by Elvis et 
al. and Ramos Almeida et al.). In addition, in the low- luminosity regime various changes in the 
structure of the AGN might occur. This is indicated by a different spectral energy distribution 
(SED) observed in these low-luminosity AGN (LLAGN)[2j, theoretically explained by a change 
in the accretion disk structure from an optically-thick and geometrically-thin accretion disk [3] 
to an optically-thin and geometrically-thick accretion flow (e.g., advection-dominated accretion 
flow; ADAF [4J) at low accretion rates (Eddington ratio r/Edd ^ 10~^). At even lower rates 
(^Edd ^ 10~^), a jet of outflowing plasma along the rotation axis of the accretion disk might 
even dominate the whole SED 0. Finally, theoretical work predicts that the dusty obscuring 
torus will disappear at low luminosities/accretion rates f6l[7]. Observational indications for its 
disappearance have been found mainly for radio-loud LLAGN [8], while remaining unproven for 
the LLAGN in general (but see also contribution by Fernandez Ontiveros et al.). 

Searches for the putative torus are most promising in the mid- infrared (MIR), where the 
thermal emission of the dust is peaking. In particular, ground-based high-spatial resolutions 
MIR observations yield valuable information about this probably compact (^ 1 pc) structure. 
It was found that the observed MIR emission of local AGN can be best explained by a torus 
consisting of many clouds [U',TI]] (and also contribution by Ramos Almeida et al.). Although MIR 
interferometry provides an excellent possibility to investigate and resolve the torus in nearby 
AGN (see contributions by Burtscher et al., Kishimoto et al., and Tristram et al.), current its 
sensitivity is not sufficient to detect most LLAGN. 

Even though lacking the high-spatial resolution, early MIR observations led to the discovery 
of a strong correlation between the nuclear MIR and X-ray emission [111 [T2] . This MIR-X-ray 
correlation can in principle be understood by a close connection between the ultraviolet and 
X-ray emission originating in the inner accretion disk region. While the hard X-ray emission 
can penetrate the obscuring torus almost unabsorbed as long as the column density remains in 
the Compton-thin regime, the ultraviolet is absorbed by the dust and re-emitted in the MIR as 
thermal radiation. Note that this explanation is independent of the torus orientation, because 
the MIR emission is nearly isotropic. Therefore hard X-rays are also a good direct tracer of 
the accretion activity and the current intrinsic power output of an AGN. By using high-spatial 
resolution MIR VLT/VISIR imaging, we could recently demonstrate for a sample of ~ 40 local 
AGN that the MIR-X-ray correlation is indeed very strong, in particular when the hard 2-10 keV 
X-ray emission is corrected for absorption [13]. The present work represents a direct extension 
of our previous work into the low-luminosity regime and is driven by the main goal of probing 
the existence of a torus-like obscuring structure in LLAGN. 

However, because of the large number of possible MIR emitters in this luminosity regime 
(Figjl]), only a first step will be presented: high-spatial resolution Gemini/ Michelle and 
VLT/VISIR MIR imaging of LLAGN. We can investigate the validity of the MIR-X-ray 
correlation for LLAGN and constrain possible contamination by circum-nuclear star formation 
affecting the MIR photometry. Star formation can significantly contribute or even dominate the 
emission at low luminosities. For a detailed description of the VLT/VISIR data we refer to |14j . 
while comparable observations are also presented in the contribution by Mason et al. 

2. Methods 

In order to separate the nuclear MIR emission in LLAGN from surrounding emission regions, 
high angular resolution data is necessary. Only this guarantees that non-nuclear contamination 
is mostly avoided, as demonstrated, e.g., in the contribution by Prieto et al. Therefore, we used 



Figure 1. Schematic 
view of the central 
structure of an LLAGN. 
The red arrows indicate 
the MIR emission of 
the various components 
and its nature. 

VLT/VISIR and Gemini/MicheUe imaging in narrow-band filters to probe the 12 /im continuum 
and polycyclic aromatic hydrocarbon (PAH) at 11.3 ^m for two samples of nearby LLAGN. These 
were selected to have absorption-corrected X-ray fluxes F2_iokeV > 10~^^erg/s/cm^, taken from 
the literature, mainly [15] and [1^. The X-ray flux limit corresponds to an estimated 12 ^umflux 
of lOmJy, assuming the MIR-X-ray correlation from [T3j, and represents roughly the lowest 
detectable flux with the ground-based MIR instruments that we used. The sample consists 
of LLAGN with a mean X-ray luminosity (logL2-iokev/ergs~^) = 40.4 at a mean redshift of 
(z) = 0.005, corresponding to a distance of (D) = 19Mpc and a resolved physical scale of 
(ro) = 40 pc. The southern sample of 18 LLAGN was observed in 2009, between April and 
September, with VLT/VISIR in the PAH2 (11.25 ± 0.59 /im) and Nellrefl (12.27 ± 0.18 /im) 
filters, and is described in more detail in [H]. These LLAGN have mixed properties, e.g., 
optical classification and are not complete in any sense. The northern sample consists of 10 
low-luminosity Seyferts drawn from the Palomar LLAGN survey [T7j and were observed with 
Gemini/ Michelle in February 2010 in the Si-5 (11.6 it 0.55 fim ) and Si-6 (12.5 it 0.6 //m ) filters. 
All observations were performed in the standard chopping and nodding mode. The data were 
reduced with the observatory delivered pipeline packages and photometric fluxes were measured 
with custom IDL routines [9]. 

3. Results 

We detected 7 out of 18 southern and all 10 northern LLAGN with all detected objects being 
visible in both filters. The high number of non-detected LLAGN can be explained by the high 
uncertainties in the X-ray luminosities that were used for the sample selection. Indeed, many 
of these have been revised to much lower values in more recent X-ray observations (further 
discussed in jl4j). Therefore we included 9 additional LLAGN with VISIR observations taken 
from the literature. All detected LLAGN appear point-like in the MIR at an angular resolution 
of ~ 0.4" , and no extended or off-nuclear emission is visible in the central 4" around these nuclei, 
similar to the findings for most of the brighter AGN [18]. The corresponding VISIR images can 
be found in p!5] and a detailed analysis of the Michelle data will be presented in a future work. 

3.1. Comparison to Spitzer - star formation contribution 

Because Spitzer /IRS spectra (angular resolution ~ 4" ) are widely used throughout the literature 
to infer AGN properties, a comparison with the Michelle and VISIR photometry is of great 




interest. In Fig. [2] we show a selection of objects from the VISIR sample. While the fluxes agree 
well for some objects, e.g., NGC 1566, the majority of objects (~ 60%) show significantly lower 
fluxes in VISIR than in Spitzer data (e.g, NGC 1097, NGC 4261, NGC 4579). In particular, 
many LLAGN do not possess PAH 11.3 /im emission on a 0.4" scale (e.g. NGC 7213). However, 
NGC 1097 and NGC 1566 most likely exhibit nuclear PAH 11.3 /^m emission as indicated by 
the comparably high PAH2 filter fluxes. We attribute the missing flux to extended diffuse 
emission regions that surround the nucleus on a 4" scale. These have a low surface brightness 
and, therefore, are over-resolved and not visible in the VISIR and Michelle images while being 
included in the Spitzer beam. In general, there is no evidence for either any correlation between 
the VISIR and Spitzer fluxes, or any trend with the optical classification. 

PAH emission is usually attributed to star formation and hence routinely used as its 
tracer. Furthermore, a tight correlation between the PAH 11.3 /im emission line flux and the 
12 /im continuum is present in local starburst galaxies [H]. Reversely we use the PAH emission 
to estimate the amount of MIR continuum emission due to star formation in the innermost 0.4" of 
the observed LLAGN. For this purpose we construct a starburst template SED by normalizing 
the IRS spectra from [19j by their PAH 11.3 /xmflux. This template SED can in turn be applied 
to the individual LLAGN by scaling with the PAH 11.3 /im emission occurring in those objects. 
With the data at hand, this PAH emission can only be constrained by using either the IRS 
spectrum or the filter containing this feature (i.e. PAH2 for VISIR), depending on which provides 
the lower upper-limit. Hence, the scaled starburst template SEDs (orange-filled curves in Fig. [2]) 
in turn represent upper limits for the star formation related MIR emission. For the majority 
of LLAGN (~ 70%), we can constrain the 12 /im continuum flux contamination through star 
formation to be less than 50% on a 0.4" scale, while for the remaining 30% the contamination 
can possibly be up to 100%. 
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Figure 2. Compar- 
ison of VLT/ VISIR pho- 
tometry (red symbols) and 
Spitzer IRS spectra (black 
line) for LLAGN. Horizontal 
error bars correspond to the 
filter pass band. Commonly 
observed emission lines are 
indicated by the dotted lines. 
In addition, the scaled star- 
burst template SED (orange- 
filled curves) correspond to 
the derived upper limit for 
the starburst contribution to 
the VISIR fiuxes. 



3.2. MIR-X-ray correlation 

As mentioned in Section 1, a strong correlation exists between the absorption-corrected hard 
X-ray and the observed 12 /im emission for local AGN. Fig. |3] shows the luminosities at these 
wavelengths for the LLAGN and the sample of local AGN from [13j. The LLAGN tightly 
follow the same correlation and extend it down to luminosities of ~ lO^^'erg/s. However, the 
observed LINERs might exhibit a minor systematic offset (~ 0.3 dex) towards higher MIR or 
lower X-ray luminosities. Nevertheless, the numbers are too low for any statistically significant 
difference, and in general, the correlation remains independent of the optical classification within 
the measurement uncertainties. A power-law fit using f itexy |20j yields: 

which is similar to previous results |13] . 
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Figure 3. Absorption- 
corrected hard X-ray luminosi- 
ties vs. the nuclear MIR lu- 
minosities for LLAGN (prelim- 
inary; filled symbols) and the 
AGN from [13j (empty symbols); 
blue squares: type 1 Seyferts 
(type 1.5 or lower); red diamonds: 
type 2 Seyferts; green triangles: 
LINERs; objects highlighted with 
central black-filled circles: "well- 
resolved" AGN from [TH] ; solid 
line: power-law fit to all LLAGN 
and AGN. 



3.3. Dependency on accretion rate 

One of the fundamental parameters determining the AGN structure is assumed to be the 
accretion rate, which can be estimated by the Eddington ratio ?7Edd = log-^boi/-^Edd, with Lboi 
being the bolometric luminosity of the AGN, .^Edd = 1-26 x 1O'^®(Mbh/^0) erg/s its Eddington 
luminosity depending on Mbh, the mass of the black hole. We estimate the black hole masses 
of the LLAGN mainly by using the Mbh-o" correlation [21] with the stellar velocity dispersions 
fj taken from the literature (mainly p2j and the Hyperleda database [23]). Furthermore, the 
bolometric luminosity is simply assumed as Lboi = 10L2-iokeV) which seems reasonable for 
LLAGN [2]. 

The Eddington ratios of all observed LLAGN and the AGN from [13] extend over a large 
range of 5 orders of magnitude, 10~^ ^ 7?Edd ^ 1- The MIR-X-ray luminosity ratio is 
constant, log(ALA(12 //m)/L2-iokev) = 0.31 it 0.36, over the whole range of r^Edd) and therefore 
independent of the accretion rate for all observed AGN. 



4. Discussion & conclusions 

The high-angular resolution MIR imaging of ~ 30 nearby LLAGN showed point-like nuclear 
emission at ~ 50 pc scale in all detected objects. No correlation to larger aperture data like 
Spitzer is evident from a comparison of VISIR photometry and IRS spectra. Therefore one has 
to be careful when deriving AGN properties from these larger aperture data. In many cases 
the observed fluxes are significantly lower and spectral features, as the PAH 11.3 /xm emission 
feature are absent on smaller scales. Furthermore, by using a starburst template SED, we 
showed that only ^ 30% of the high-resolution data may still be dominated by MIR emission 
from star formation occurring on ^ 50 pc scale. The majority {'^ 70%) of LLAGN are, however, 
dominated by AGN emission on this scale, while its exact origin - torus, jet or accretion disk 
- remains undetermined at this point. Furthermore, the MIR-X-ray correlation is valid for all 
AGN observed at high-angular resolution in the MIR over the whole probed luminosity range 
from 10^*^ to ~ 10^^ erg/s. Owing to the fact that no dependence of either optical classification, 
or radio-loudness, or host morphology was found, this correlation might put strong constraints 
on the universal AGN structure. Interestingly, none of the theoretically predicted changes in the 
accretion structure with decreasing accretion rate, e.g., onset of ADAF or the jet-dominance, 
seem to have any imprint on the MIR-X-ray luminosity ratio. However, the sample is likely 
incomplete (selection based on availability of 2-10 keV data) and possibly biased against highly 
obscured (Compton-thick) objects. This along with the small total sample size prevents us 
from reaching any final conclusions on the AGN structure and nature of the MIR emitter in 
the low-luminosity regime. At this point we can only conclude that the results agree with 
those for brighter AGN, for which the MIR emission is well explained by a dusty clumpy torus. 
Therefore, such a torus might also exist in LLAGN, i.e. the standard unification model of AGN 
might still hold in this regime, in agreement with [16^124). Furthermore, these results imply that 
in LINERs, the emission is AGN-dominated in at least central tens of parsecs (compare |25|). 
Finally, the MIR-X-ray correlation represents a powerful observational tool to convert between 
both wavelength bands for any AGN. In order to settle the question about the nature of the 
MIR emission, SED modeling with particular focus on the high-resolution infrared data seems 
most promising (see also contributions from Prieto et al. and Fernandez Ontiveros et al.). 
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